Scaling and Interference in the Dissociation of Halo Nuclei 
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The dissociation of halo nuclei through their collision with light and heavy targets is considered 
within the Continuum Discretized Coupled Channels theory. We study the one-proton halo nucleus 
*B and the one-neutron halo nucleus ^^Be, as well as the more normal ^Be. The procedure previously 
employed to extract the Coulomb dissociation cross section by subtracting the nuclear one is critically 
assessed, and the scaling law usually assumed for the target mass dependence of the nuclear breakup 
cross section is also tested. It is found that the nuclear breakup cross section for these very loosely 
bound nuclei does indeed behave as a-\-bA^^^. However, it does not have the geometrically inspired 
form of a circular ring which seems to be the case for normal nuclei such as ^Be. We find further that 
we cannot ignore Coulomb-nuclear interference effects, which may be constructive or destructive in 
nature, and so the errors in previously extracted B(E1) using the subtraction procedure are almost 
certainly underestimated. 

PACS numbers: 25.60.Dz, 25.70.Dc, 24.10.Eq 



The study of electromagnetic dissociation of halo nu- 
clei is an important area that supplies invaluable infor- 
mation about their multipole responses and consequently 
the details of their structure 0, Q . The theory usually 
employed for the purpose is the relativistic Coulomb ex- 
citation theory developed by Winther and Alder P, Q 
which hinges on the Fermi- Weisacker- Williams (FWW) 
virtual photon method. Since the measurement of the 
dissociation cross section supplies the combined Coulomb 
and nuclear contributions, one is forced to subtract the 
latter. The common prescription employed for this sub- 
traction procedure is the so-called scaling law: the nu- 
clear breakup cross section should scale linearly with the 
radius of the target, and thus by measuring the cross 
section for the breakup of the halo nucleus in the pre- 
dominantly nuclear field of a light target, one attempts 
to extrapolate to heavy targets assuming the validity of 
the scaling law. The Coulomb dissociation cross section 
for the halo nucleus on the heavy target is then simply ob- 
tained by subtracting from the experimental cross section 
an extrapolated nuclear one calculated according to some 
prescription. It is this "nuclear-free" cross section which 
is fitted by the FWW result in order to extract the B{EL) 
distribution |llllE[l^l,iail,i2„i3,i4 15,16,13]. 

Two factors convinced us to criticall y a ssess this pro- 
cedure. Several theorists [11 [H H M El 113 have 
recently cast doubt on the relative importance of the nu- 
clear contribution to the dissociation cross section, claim- 
ing in some cases that this contribution can be signifi- 
cantly larger than the Coulomb contribution in as heavy 
a target as lead, invalidating the scaling law. The sec- 
ond important factor is the need to supply a quantitative 
assessment of the Coulomb-nuclear interference terms in 



the cross section. 

The purpose of this letter is to settle these issues by 
performing a full Continuum Discretized Coupled Chan- 
nels calculation for *B, ^^Be and ^Be dissociations in the 
fields of light-, medium- and heavy-mass targets at three 
laboratory energies where data are available. 

Coupled channels calculations were therefore per- 
formed to calculate the elastic breakup (also called 
diffraction dissociation) arising in a three-body model 
consisting of a two-body projectile incident on an inert 
target [24l [25l | . Between the target and each of the two 
projectile components we specify optical potentials whose 
imaginary parts describe the loss of flux to channels be- 
yond elastic breakup. We calculate all orders of the tidal 
effects of the optical potentials as they deform and break 
up the projectile. 

For both *B and ^^Be, projectile states were included 
for relative motion up to partial waves £max — 3 and 
energy emax = 10 MeV. This energy range was divided 
into 20 bins when I — Q for ^B and when £ = 0,1 for 
"Be. There were 10 bins when ^ = 1, 2 for ^B and ^ = 2 
for ^^Be, and 5 bins in the remaining partial waves, all 
evenly spaced in momentum k. The coupled channels 
for the scattering of the projectile on the target were 
solved up to i?max = 500 fm, and for partial waves up 
to if -Rmax where K is the wave number for the incident 
beam. The optical potentials have the same parameters 
for all energies, as given in Tabled using radii calculated 
with an A^^^ contribution from the target mass number. 
We thus impose a regular target scaling in the radial 
geometry of the component-target potentials, in order to 
examine the contributions from varying target size and 
varying dynamical conditions. 
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TABLE I: Optical and binding potential parameters for *B 
and ^^Be projectiles. 
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The results of the integrated nuclear breakup cross sec- 
tions obtained from the CDCC calculations are shown in 
figure n circles and squares for ®B at Eiat = 44 and 
70 MeV/n respectively, left-triangle, down-triangle and 
right triangle for "Be at Eiab = 44, 70, and 200 MeV/n 
respectively, and full diamonds for ^Be at 100 MeV/n 
(after [i^). These are plotted as functions of the cubic 
root of the mass of the target nucleus, proportional to the 
target's radius. They are all rather accurately accounted 

1 /3 

by a linear dependence on Aj( , but with varying coef- 
ficients. We might expect that the total direct reaction 
cross-section is proportional to the area of a circular ring 
given l27| by: 



an « 2Tra{Rp + Rt) 



(1) 



where a is approximately the diffuseness of the optical 
potential that describes the elastic scattering, and Rp + 
Rt is the sum of the projectile and target radii. In |27l |. 
it was found the Eq. ^ describes the data quite well. 
Guided by Eq. iQJ we anticipate the following form of 
the nuclear breakup cross-section: 



CTAT = Pi + P2A. 
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(2) 



where the parameters Pi and P2 depend on the projectile, 
and may also depend on the bombarding energy and the 
structure of the target. We performed linear fits to the 
nuclear breakup cross sections as a function of A^^^ to 
obtain the values for Pi and P2 shown in the legends of 
Figs. U] 

The nuclear breakup cross section calculated with 
CDCC for ®B and ^^Be do show the A^p^^ dependence 
of Eq. ||2Jl as seen in figures ^ but do not always follow 
the form given by Eq. in particular as, in most cases, 
the cross section fits have Pi < 0. At this point we should 
mention that our results disagree completely with those 
of |23| , where the nuclear breakup cross section was found 
to scale as At- (In fact such a behaviour could only be 
possibly true for a very weakly interacting system where 
the whole volume of the target could be effective.) 

This is not the case for ^Be, a normal non-halo nucleus, 
where scaling holds and Eq. ^ is fully satisfied, where 
both Pi and P2 are positive as seen in figure ^ The 
different behaviour in the scaling of the nuclear breakup 
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FIG. 1: Elastic nuclear breakup cross-section for ^B, ^^Be 
and ^Be projectiles at the indicated energies, as a function of 
target mass number At, along with linear fits. 

cross-sections for ^B and "'^^Be could well be a manifesta- 
tion of their halo nature. 

The fact that the nuclear cross sections do not exactly 
satisfy Eq. with Pi proportional to the radius of the 
projectile leads us to conclude that nuclear breakup must 
be estimated from realistic dynamical models, not from 
Eq. QJ. As well as the present CDCC model, Glauber 
few-body procedures can be used, as done in refs. [23.12^ 

m. 

The analysis of experimental data [2^ |2^ [s^l started 
with the following expression for the breakup cross sec- 
tion: 



da 

dE* 



= S 



.doc 
dE* 



L{At) 



daj^^C) 
dE* 



(3) 



where S is the ground state spectroscopic factor, E* is 
the excitation energy. The dac /dE* is from the Coulomb 
FWW virtual photon formula, and the da{^^C)/dE* is 
what is observed for a ^^C target. The L{At) is a scaling 
factor, which may be determined either by fitting Eq. © 
to the data, exploiting the different shapes of the excita- 
tion spectra, or by calculating 



L(^T)-<(AT)/aJv"(12) 



th/ 



(4) 



using eikonal calculations of nuclear breakup a)^{AT)- 
The ^^C is simply a reference nucleus. Of course 
the above incoherent sum ignores completely Coulomb- 
nuclear interference effects, which we discuss further be- 
low. 

In m El the reaction "Be + ^ospb at E = 520 
MeV/n is considered. By adjusting S and L, these au- 
thors obtain P(208) — 5.6 ±0.4. This is close to our value 
of L, which can be extracted from Fig. □ at S = 200 
MeV/n, namely L = 5.9. Our results also agree with the 
eikonal calculations of [sj. In fact, a.t E = 200 MeV/n, 
it has been further confirmed to us [33 that within the 
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FIG. 2: Total breakup, Coulomb only, and nuclear only con- 
tributions for ^^Be projectile as a function of A)^^ , for a lower 
radial cutoff b^nm ~ 20 fm. 
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FIG. 4: Total breakup. Coulomb only and nuclear only 
breakup cross-section for ^Be projectile 
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FIG. 3: Total breakup. Coulomb only cross-section and nu- 
clear only breakup cross-section for *B projectile, for a lower 
radial cutoff femin = 20 fm. 



eikonal formalism of , the value of L comes out about 
5.5, very close to our model result. 

On the other hand, using the same incoherent cross- 
section formulae, Ref. 3Q| analyses the reaction ^^Be -|- 
208pb at £; = 70 MeV. These authors found for L the 
value 2.1 ± 0.5, much smaller than our result L = 3.82 
(Fig. ^ but in accordance with the geometrical value 
obtained from Eq. (^3), which we have shown to be not 
valid for halo nuclei for at least some halo nuclei. It is 
conceivable that the reason for this discrepancy resides 
in the neglect of the interference terms, which, as we see 
next, can be rather large. 

The Coulomb and nuclear potentials combine together 
coherently to give breakup, and the destructive inter- 
ference between these potentials in the surface region 
is well known. If, however, the nuclear and Coulomb 
breakup cross sections contribute largely to different par- 
tial waves, then the total breakup cross sections will add 
incoherently and Eq. |j2Jl should be accurate. 

In order to answer this question definitively at least in 
our test cases, we have performed further CDCC breakup 
calculations with both Coulomb and nuclear transition 
potentials as in [ssl 0| , along with sufficient radial and 
partial wave limits to encompass all the resulting breakup 
cross sections. For reference we also performed 'Coulomb 



only' calculations, where there are no nuclear potentials 
at all. Any pure Coulomb calculation needs at least to be 
limited by a minimum impact parameter 6min to simulate 
the effects of nuclear absorption during grazing and closer 
collisions. 

The combined Coulomb and breakup calculation will 
therefore give a total 

o'cw = + ctat + (J/ (5) 

which defines an interference term cr/ by the difference 
with the sum of Coulomb-only and nuclear-only calcula- 
tions. We have found that uj is sometimes negative (de- 
structive interference), sometimes positive (constructive 
interference), and often large. Thus, although one can 
construct ctat from some scaling model, the mere sub- 
traction of it from the data would give a "contaminated" 
Coulomb breakup cross-section: 

O'C = CTc + CT/ = (JcN — ctn ■ (6) 

The use of, say, the equivalent photon method as done 
in 0i 'E'l to extract the B{El) distribution for and 
^^Be or for that matter any other halo nucleus could be 
questionable if ai is large. To ascertain typical sizes of 
tj/, we first show in figures [3 El and ^ the cross-sections, 
fCj fiv and ucn as the short dashed, long dashed and 
solid lines respectively, for the ^^Be and ^B at Eiab — 44 
MeV/n, and ^Be at Eiab = 100 MeV/n, using 6„iin = 20 
fm in each case. 

To see this, in Figs. [21 and O we plot as the dot-dashed 
line the values of gcn ~ ctat, which should 'ideally' agree 
with ac as the short dashed line. We also plot as the 
dotted line the incoherent sum uc + iyN, whose difference 
from the ac n solid line indicates the effects of interfer- 
ence. It is evident that, for large At, the interference 
term is destructive for *B (the solid line acN < ctc+n 
dotted line), so its neglect may lead to unrealistically 
smaller B{E\) distributions. Conversely, it is construc- 
tive for "'^^Be, yielding an unrealistically larger B{E1). 

Experimentalists (e.g. 17]) often try to minimise these 
interference problems by restrictions to small excitation 
energies, and/or to large impact parameters. We focus 
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FIG. 5: Ratios of the true to the contaminated Coulomb 
breakup cross sections ac/o-c = o-c/icrcN — crjv) as a func- 
tion of the lower radial cutoff bmin, for four different targets. 
Results for *B are shown in the upper panel, and for ^^Be in 
the lower panel. 

on the impact parameter restrictions, implemented by 
means of a maximum angle 0,nax for integrating the cross 
sections over the centre of mass angle of the projectile 
fragments. By semiclassical theory, this angle is related 
to fcinin by f?max = 2?7/ (fc^min). We therefore use the same 
restrictions on the calculated breakup cross sections. In 
Fig. jSjwe show the ratios of the 'true' to 'contaminated' 
Coulomb breakup cross sections acjoc = od ipcN — 
CTjv) as a function of the lower radial cutoff &min- 

Ideally the ratios should be unity, but in fact we see 
that rarely do these values even tend to unity for large 
^min- Only for the projectile on ^°^Pb does this oc- 
cur, and the deviations from unity are worse in the ^^Be 
case. These deviations from unity indicate either that the 



long tail of the ^^Be ground state wave function gives 
rise to small but significant deviations from the pure 
Coulomb results even at impact parameters ^ 30 fm, 
or that diffraction effects are large enough to break the 
semiclassical connection between &min and ^max- Equiv- 
alent plots for the ^Be projectile (not shown) give ratios 
far from unity, from diffraction or refraction causing nu- 
clear breakup fragments to come out at very forward an- 
gles, implying that there is no 'safe' angular region where 
Coulomb effects dominate |2^. In general, it is clear that 
theoretical breakup calculations of Coulomb-nuclear in- 
terference are needed for accurate results for the breakup 
of ^^Be and ^Be on any target. 

In conclusion, we have given evidence through detailed 
CDCC calculations of the scaling behaviour of the nu- 
clear breakup cross-section. This cross-section does scale 

1 /3 

as A^'", the mass number of the target nucleus, but the 
scaling does not follow the geometrical form as in normal 
nuclei. The nuclear contribution can be as small as 1/15 
of the values calculated in refs. 0,0,|23|. We have fur- 
ther calculated the "error" due to the nuclear-Coulomb 
interference in the extracted B{E1) distribution if the 
subtraction acN—o'N is employed in conjunction with the 
virtual photon method. We believe that a full quantum 
calculation, including both Coulomb and nuclear poten- 
tials on equal footing, such as CDCC, is required to get 
credible numbers for the B{E1) distribution of ncutron- 
and proton-rich (as well as some stable) nuclei. 
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